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The purpose of this research is to investigate specific 
processes that contribute to large scale equatorial SST 
variability, and by comparison with observation, verify the 
realism of the model's thermal response to prescribed 
forcing in the atmosphere. Building on previous work, more 
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ocean part of Rennick's (1985) coupled ocean-atmosphere 
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forcing is conducted. 


The dynamic ocean model is based on the shallow water 
momentum equations forced by a surface wind stress. Lt 35 
linearized about a motionless state with a zonally sloping 
pycnocline depth which is in balance with the surface stress 
caused by the zonal wind of the atmospheric basic state. 
This study investigates the effect of turbulent vertical 
mixing of heat, in contrast to horizontal advection , for 
the generation of SST anomalies in the equatorial region of 
the Pacific. Results showed that the SST anomaly produced 
by turbulent mixing was two orders of magnitude smaller 
than, and 90 degrees out of phase with, the SST anomaly 


generated by horizontal advection. 


Although the small magnitude of the anomalous entrainment 
heat flux would seem to justify its neglect, the phase 
difference raises the question of whether it would be 
significant in a coupled ocean-atmosphere system. 
Therefore, two cases of the coupled model were designed to 
investigate if the different SST responses caused by turbu- 
lent mixing and by horizontal advection would cause growth 
of a coupled disturbance. The results were that neither 


ocean process causes growth. 
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I. INTRODUCTION 


A. BACKGROUND 


The El Nino/Southern Oscillation (ENSO) is the dominant 
global climate signal on time scales of a few months to a 
few years (Rassmusson, 1984). ENSO's have appeared every 
few years for centuries, and there have been at least six 
since 1951. 

The aspect of an ENSO event is that once it starts, it 
evolves rather predictably for as long as 18 months. The 
meteorlogy of a typical ENSO event can be described in terms 
of the spatial evolution and the time of occurrence relative 
to the seasonal cycle, which are strikingly similar for most 
ENSO events (Rassmusson, 1984). The year prior to the event 
is characterized by a strong trough over the Maritime 
Continent (Australian, Indonesian, Southern Asia Monsoon 
regions) and a strong south Pacific ocean High. This pres- 
Sure gradient is called the Southern Oscillation Index (SOT) 
and is defined as the difference in surface pressure between 
Tahiti and Darwin. A slackening of the pressure gradient 
(SOL) at the end of the first year, along with a weakening 
of the easterlies, causes positive sea-surface temperature 
anomalies (SSTA'S) to propogate eastward across the equato- 
rial Pacific Ocean to the coast of Chile by one or more of 


the following possible dynamic mechanisms (Philander, 
OSD): 


(1) anomalous southward flow of warm water due to the 


slackening of southeast trades; 


(2) changes in large-scale equatorial zonal currents 
associated with local changes inthe curl of the 


Wind stress; and 


(3) local and remote warming due to weakening of equato- 

rial easterlies. 

The ENSO year is characterized by rapid positive and 
negative surface pressure changes at Darwin and Tahiti and 
sea-surface temperature increases near the Peru coast. Also 
of crucial importance is the major eastward shift of the 
Maritime Continent convection region. In the upper tropo- 
sphere, an anticyclonic couplet straddles the region of 
enhanced rainfall in the Central Equatorial Pacific, causing 
a stronger-than-normal eastern North Pacific jet wa 
persistent upper-troposphere easterlies. At the onset of El 
Nino, the Southern Hemisphere experiences anomalies in large 
scale flow in the formation of a huge anticyclone at 15 
degrees south. Northward and eastward of the anticyclonic 
center is an apparent wave train of positive and negative 
anomaly centers which lies along a great circle and may 
reflect a Rossby wave train (Wright, 1977). In the western 
North Pacific ocean, low-level westerlies extend eastward 
during the ENSO event and create an environment of low level 
cyclonic relative vorticity near the dateline. 

Scientists thought they had a reliable picture of how an 
ENSO event started and developed, but then the 1982-1983 


ENSO event occurred. Even though in scientific terms it 
came along at the expected time, it did not occur in EiS 
expected manner. The warm water typical of ENSO usually 


first appears off the coast of South America early in the 
year and then gradually spreads westward. Instead, the 
warming in 1982 began in May in the central Pacific and 
gradually spread eastward. By strict definition this was 
not an ENSO event, but it turned out to be the strongest 
ENSO event in nearly a century. 

The 1982-1983 ENSO event is being intensely studied, 
because it was the most intensely measured event. It was 


found that the sea-surface temperature (SST) vice SSTA were 


10 





closely related to the anomalous convective activity. The 
eastward migration of the region of maximum SST was accompa- 
nied by an abnormally large heat loss to the atmosphere, and 
a westerly wind anomaly over and to the west of the SST, 
maximum. The combined system of warm SST, convection and 
westerly wind anomalies then progressed steadily to the 
east. This ENSO event had two sequential but interrelated 
events (Rassmusson, 1984). The first event was the initial 
warming near the eastern boundary associated with the weak- 
ening of the eastern South Pacific ocean anticyclone, and 
the second was the broad scale warming of the equatorial 


Pacific from near the date line eastward. 


ENW THE PROBLEM 


Numerical models based on the shallow water eguations 
have been very useful in investigating and studying ENSO 
events. Hurlburt et al. (1976) and McCreary (1976) inves- 
tigated ENSO-like events by forcing their simple models with 
abrupt reductions in surface wind stress. Busalacchi and 
O'Brien (1981) and Busalacchi et al. (1983) successfully 
simulated the observed variability of sea level for a period 
1961-1978 using a shallow water equation ocean model with a 
prescribed surface wind stress based on monthly mean aver- 
ages. The response of the tropical atmosphere to anomalous 
equatorial heat sources using the shallow water equations 
was studied by Gill, (1980). Gill and Rasmusson (1983) used 
a heat source based on the observed mean of outgoing long 
wave radiation from June 1982 to January 1983 and found 
there was an agreement between their model wind and the 
corresponding observed wind field anomalies. 

The degree of interaction between the simultaneously 
evolving atmospheric and oceanic elements of ENSO events has 


been explored in a limited sense by Lau (1981) and by 


]EP 


McCreary (1983). In their models the atmosphere assumes one 
of two prescribed states dependent upon the distribution of 
the model-determined pycnocline depth anomaly. The coupled 
model produced an oscillation between ENSO (warm ocean) and 
a normal (cold) ocean. 

The coupled air-ocean models were designed to answer 
critical questions concerning the initiation of an ENSO 
event. The questions involve the mechanisms that are 
responsible for ENSO, the timing of the event relative to 
the annual cycle, the time interval between events, and the 
eastward migration of the event. Rennick (1985) utilized a 
simple numerical model which allowed interactive coupling 
between the atmosphere and the equatorial ocean. An inves- 
tigatlon was made to determine the extent to which eastward 
extension of the warm water pool in the western Pic iiS 
during ENSO could be attributed to the atmospheric circula- 
tion induced by the warm water itself. These models ignored 
vertical mixing and attempted to isolate the role of hori- 
zontal advection of temperature. The result of this numer- 
ical experimentation was that the atmospheric circulation 
due to an anomalously warm western Pacific did not induce a 
self perpetuating eastward migration of warm water and 
anomalous convection. 

The present problem involves using the same dynamic 
models used by other investigators (ie. McCreary, l9 
1983; Busalacchi and O'Brien, 1981; Bussalacchi et al., 
1983; Lau, 1981; Rennick, 1985). Rennick and Haney (1985), 
using a linear coupled atmosphere-ocean model, found ERIE 
the analytic solutions are very sensitive to whether heating 
in the atmosphere is parameterized in terms of the oceanic 
pycnocline depth anomaly (Philander et al.,1983) or in terms 
of the advectively determined sea surface-temperature 
anomaly (Rennick and Haney 1985,Rennick 1985). The objec- 


tive of this research is to investigate what determines sea 
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surface-temperature variability using the ocean model part 


of the coupled atmosphere and 
Rennick (1985). This will be done by incorporating more 


equatorial ocean model of 


ocean thermodynamics than previously considered. 
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II. IHEORTES MODELS 


A. INITIAL COUPLED AIR-SEA MODEL USED BY RENNICK (1985) 


The numerical model used here is a coupled model in 
which the atmosphere and ocean are each represented by the 
shallow water equations. The model will be used primarily 
in the ocean-only mode, and, therefore, all theory pertinent 
to the ocean will be discussed in detail. 

The equations for the ocean are solved in an equatorial 
basin having 120 degrees longitude by 30 degrees latitude 
centered on the equator at 140 degrees west longitude. The 
ocean is forced by a surface wind stress where the mean wind 
is easterly at five meters per second, and an anomalous 
westerly wind patch (Guassian) is utilized to depict the 
anomalous conditions experienced during ENSO. 

Two important features are included in this model which 
are absent in earlier work. First, the annually averaged 
pycnocline depth is set at two hundred meters in the west 
and eighty meters in the east. This zonally sloping pycno- 
cline depth is in balance with the surface stress attribu- 
table to the mean zonal wind. This modification allows the 
pycnocline anomalies computed by the model, especially in 
the eastern part of the basin, to be smaller than those 
associated with a uniform mean pycnocline depth. 

Secondly, rudimentary mixed layer physics is included 


along with the dynamic processes of the shallow water equa- 


tions. This allows the the sea-surface temperature to be 
computed directly. This is very important in an interactive 
system, as pointed out by Hanson (1983), and Rennick and 


Haney (1985). 


The model equations for the ocean are then written as 
(1) du/dt - fv * g' dh/dx - Fx * KV*u 
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Cidade + bs E dh/dy = Fy + K Váwv 


(3) dh/dt + NY* (Hv) = O 


(â) dT/dt « N7* T> -X 


H(x) is the pycnocline depth of the motionless ocean basic 
state, T is the mean sea-surface temperature, and g is the 
reduced gravity for the two layer ocean. There is no mean 
flow in the ocean. Fx and Fy are the respective wind stress 
in the eastward and northward directions. The ocean thermo- 
dynamic equation (4) includes horizontal advection by the 
dynamically produced currents and a surface heat loss term, 
XT. Based upon arguments of Schopf (1983), the time scale 
for the surface heat loss (1/X) is chosen to be about one 
year. In this initial model, effects of vertical mixing in 
the ocean were ignored in order to isolate the role of hori- 
zontal advection. 

The model's basic state is defined by solid (constant 
angular velocity) eastern rotation in the atmosphere 
(U=-5m/s), the slanting pycnocline, and the mean sea-surface 
temperature T decreasing linearly from 30° C at the western 
ocean boundary to 25° C at the eastern boundary. In order 
to define the initial (anomalous) conditions for the ocean 
model, a sea-surface temperature anomaly was first imposed 
at the western edge of the ocean. This sea-surface tempera- 
ture anomaly was held fixed while the atmosphere was spun up 
from its basic state (zero anomaly) to one of approximate 
equilibrium with the heating anomaly associated with the 
imposed sea-surface temperature distribution. The resulting 
zonal wind anomaly was used to compute an approximately 
balanced mixed layer depth anomaly. This mixed layer depth 
anomaly together with the original SST anomaly constitute 
the oceanic initial conditions. The currents were initially 


Zero. The wind forcing was the total (solid rotation plus 


1) 


anomaly) wind after spinup. In the ocean-only mode the 
atmospheric surface wind field is held fixed, and the ocean 
is allowed to evolve from its initial state in response to 
the fixed surface stress. Since the surface wind field was 
held fixed, there is no feedback due to changes in sea 
surface-temperature. 

The ocean's response to these conditions (Rennick, 19688 
was an eastward propagating Kelvin wave whose currents 
produced warm advection along most of the equator. When the 
Kelvin wave reached the eastern boundary (80 degrees west 
longitude), some of its energy reflected westward as a 
Rossby wave. As this Rossby wave migrated westward, it 
caused the sea-surface temperature tendency to reverse, 
thereby defining the peak of the warm event in the model 
ocean. Most important, ,however, is the fact that the 
resulting SST anomaly was guite weak. These results point 
out the need for further study, including the possible need 
to include ocean to atmosphere feedbacks, and additional 
physical processes such as vertical mixing and radiation in 
the determination of the sea-surface temperature. These 


results provided motivation for the present thesis. 


B. INCLUSION OF WIND STRESS ANOMALY FORCING 


The additional oceanic physical processes were tested in 
the ocean-only version of the model forced by a prescribed, 
but moving, wind stress anomaly field characteristic of an 
ENSO event. The wind patch used to better describe the ENSO 
wind field was suggested by Blundell and Gill (1983). 

A simple analytical form was chosen to model the 
observed stress data. The form of the eguation denoting the 


eastward component of stress (Fx) is 


(5) Fx = F(t) e G(KST) ligo 
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F(t) is the time envelope for a forcing of constant 
shape, propagating with constant speed in the X-direction. 
The wind stress(Figures 4.1-4.6) lasts for fifteen months, 
with a sinusoidal rise and fall. 

G(KSI) is the wind forcing patch and covers sixty 
degrees of longitude. CEST) is given by 


"H 


(6) G(KSI) cos(3 * KSI) for |KSI| < 77/6 


(7) G(KSI) = 0 for |KSI| > 7776 


where KSI is expressed in radians for the computer model and 


is given by 
GM KSI = X = CE. 


Here X is longitude (radians) measured eastward from the 
western boundary, t is time and c is chosen to be a motion 
of 90 degrees longitude in ll months (Rennick,1985). 

The north/south extent of the wind patch is described by 


the function H(y), where y is latitude (radians), and 
(9) H(y) = (1 - theta”) exp (-3 + theta*/4) 
(10) theta = ay 


The value of (a) is 5.73, in order for H(y) to have a zero 
crossing at 10 degrees latitude. In the model, if H(y) in 
equation (9) is less than zero, then H(y) is set equal to 
Zero. This was done because values dropped off very 
quickly, and an easterly wind forcing over any part of the 
ocean domain was undesirable. 

Figures 4.1-4.5 show a comparison of the observed anoma- 
lous wind fields for T3500 mb from mn Oil OS EA 
Nino/Southern Oscillation Quick Look Atlas (1983) with those 
of the computer generated anomalous wind fields at compa- 
rable times (lower panel). The lower panel encompasses an 


area from 15º north to 15° south latitude and 130° east to 


EZ 


50° west longitude. The central longitude is at 140° west 
vice 180° longitude (upper panels ) and therefore a box has 
been included in the upper panels showing the basin bound- 
aries of the model. The magnitude of the wind fields differ 
because of the parameters selected for the model, but the 
direction of the winds and the movement of the patch corre- 
Spond very well with the observed field. Figure 4.6 is a 
time-longitude plot of the computer generated anomalous 
surface zonal wind at the equator and shows its eastward 
movement, magnitude and location for five hundred days. 
Divergence and convergence patterns can be inferred from 
Figure 4.6 in that convergence will occur ahead of the wind 


patch and divergence will occur behind. 


C. INCLUSION OF OCEAN MIXING MODEL 


To investigate specific processes that contribute to 
large scale eguatorial SST variability, an integration of 
the temperature eguation is integrated over depth. Assuming 
that non-linear terms are negligible, the eguation has the 


form: 


(11) dT'/dt = -V *VT' - V' "VT - H'/H (WT), - (WT) 
) - (WT)'/H * (WT)' /H ü 


The overbar represents a mean value and the prime represents 
an anomaly. The term (WT)' represents the turbulent 
vertical heat flux, with (WT) the climatological value. The 
subscript O refers to the sea-surface and -H to the base of 
the mixed layer. The first term is the change of anomalous 
temperature with time. The second term represents the hori- 
zontal advection of the SST anomaly by the mean current and 
on the basis of available data (Niiler, 1982; Weare, et.al 
1981; Wyrtki, 1981) can be considered negligible. The third 


term is the horizontal advection of the mean SST by the 
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anomalous current. This" the term Gill (1983) and 
Harrison and Schopf (1984) found to dominate the SST anomaly 
development during the early phases of ENSO. Term three is 
the only term which has been utilized in Rennick's (1985) 
studies. 

Present estimates suggest that some of the other terms 
may become Significant as a warm anomaly moves eastward into 
a region where the pycnocline depth is smaller. One of 
these terms is term four of equation (11). This is a SST 
anomaly generating term, where anomalous warming or cooling 
occurs by converging the normal amount of heat into an 
anomalously thin (H'<0) or thick (H'>0) layer. Term five of 
equation (11) is the anomalous surface heat flux and 
includes all components of heating such as sensible, latent, 
and solar heating. It is term six, which is the anomalous 
entrainment heat flux, that will be of most interest in this 
thesis. 

Conducting a rough estimate of the scale of all the 
terms, it can be shown that no term is larger than the 
anomalous entrainment heat flux term, except possibly the 
term which is already included in the model. The horizontal 
advection of mean SST by the anomalous current is the most 
important, but it is believed that the contribution of the 
anomalous entrainment heat flux is perhaps comparable. 


Therefore, equation (11) will be simplified as follows: 
(12) dT'/dt - -U'* dT/dx « (WT)' /H 


This equation states that the change of anomalous tempera- 
ture with time is equal to the horizontal advection of the 
mean SST by an anomalous east-west current plus the anoma- 
lous .entrainment heat flux divided by the mean pycnocline 
depth. 

To close the ocean model it is necessary to include a 


consistent mass flux when entrainment heat flux is present. 


ILS, 


Integrating the continuity equation for an incompressible 


fluid gives 
(13) dH'/dt = -Y° (VH) = We' 


H' is the anomalous pycnocline depth, H is the mean 


— A 
pycnocline depth, and HV is expressed as 
O 
— A 
(14) HV = IE 
-H 


We' is the anomalous entrainment mixing and is expressed 


in terms of the entrainment heat flux as 
(15) We' = - (WT) | PZA M 


AT is the temperature jump at the base of the layer and 
is prescribed as 3* C in the model (Wyrtki, 1981). 

In order to obtain the anomalous heat SINE at the base 
of the layer a modified version of the Kraus and Turner 
(1967) bulk closure model of the upper ocean will be used. 
In this model, 


(16) (WT), = -As ((((1+K) * m * (WSTAR**3))/ (ALPHA * 
g * H)) + (WT), - (1-K) *As ((WT), } - S(F,H)} 


The left hand side is the total (mean plus anomalous) 
entrainment heat flux. It is necessary to subtract the mean 
entrainment heat flux to get the anomalous entrainment heat 
fluk: "AE the heavy side function. WSTAR is the friction 
velocity in the water due to the surface wind stress, ALPHA 
is the thermal expansion coefficient for sea water, and g 1s 
gravity. 

The coefficient K—0.2 (Stage and Businger, 19810 
accounts for dissipation of turbulence produced by shear and 
convection in the upper ocean, while the coefficient m—2.0 
(Niiler and Kraus, 1977) accounts for the fact that the 
production of turbulence by stress driven shear in the upper 
ocean is only proportional to(WSTAR **3)/ (ALPHA * g œ H) 
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(see Niiler and Kraus, 1977). S(F,H) is the stabilizing 


effect of solar radiation where E 


(17) S(F,H) 2 2 + (.5 + (F,- E, ) - (1./H) . Jr(z)dz) 


-H 
F(z) is the downward flux of solar radiation. ME JES 
exponential in form, and it is a maximum at the surface and 
minimum at greater depth. Since F(z) is not linear in Z, 


but rather exponential, its effect is to stably stratify the 
mixed layer, and this requires energy from the wind and 
convective mixing in order to mix the temperature to an 
isothermal state. 

Equations (16) and (17) are used to obtain the solution 
for the anomalous heat flux. The anomalous form is derived 
under several assumptions. The climatological state is 
assumed such that (WT) „>O and (WT), «0. Also,in the general 
state (WT )¿>0. Thus, the turbulent flux at the surface is 
always upward and the mean turbulent flux at the base of the 
layer is downward. These assumptions cause A ((WT) 35(WT),, 


so equation (16) becomes 


(18) (WT), = -A* ((((1+K) * m * (WSTAR**3))/(ALPHA * 
em ee (WI — S(F,H)) 


The next step is to write every term in equation (18) 


into a mean plus an anomaly, 


(19) (WD + QT) EEK) *_m * 
KE TOP se ea Ke (WI). = S(F,H) 8 
(((1+K) © m + (WSTAR**3/H)')/(ALPHA + g )) +K + 
Cy (F 


and then substitute into equation (19) for the mean terms in 
the brackets on the right hand side and subtract wr. from 
both sides to get the following form for the anomalous 


entrainment heat flux. 


Zu 


(200 (wr), = -A* COD OE E 
(WSTAR**3/H)')/ (ALPHA " g ) + K» (WT), - «Sl 
(F,H)) - (WT), | 


Expanding the anomaly of a ratio! as in the 
(WSTAR**3/H)' term in equation (20), and neglecting the 
other terms (surface heat flux and solar radiation), equa- 


tion (20) can be written as 


(21) (WT), = ÃO CT * (IE) + mos 
(WSTAR**3))/ (ALPHA . g a H )) > 
(((WSTAR**3)'/(WSTAR**3)) - (H'/H))) - (WT) 


The final form necessary for the computer mixing model 
nescessitated one more assumption. The climatological 
entrainment heat flux at the base of the layer balances much 
of the surface warming and is of the following form (Wyrtki, 
1981), 


(22) (WT), = - .75 * 00, 


where QO is the total downward heat flux at the sea-surface 
under climatological conditions. This means that (WT) is 
generally negative, causing active entrainment and a down- 
ward flux of heat at the base of the layer in a long term 
climatological sense. 00 is the net surface heating of the 
ocean and is obtained from Weare's Atlas (1980). Using a 
least squares fit on the mean climatological data, 00 was 
found to be 48.14 W/m? at the western boundary and 83.00 
W/m? at the eastern boundary. The computer model required 
QO(x) for 201 grid points along the equator, and therefore a 
straight line was fitted between the above two boundary 
values. 

Substituting? (22) into (20), gives the final form of the 


anomalous entrainment heat flux as 
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SD AN 79. = 00... (((I*K). ^ * m o 
ABS(U**3))/(ALPHA * g © H)) * (((RHOA + 
CD)/RHOSEA)**1.5) * (((3 * U')/U) - (H'/H))) + 275 
è QO 


RHOA is the air density, RHOSEA is the density of sea water, 
and CD is the drag coefficient. U' is the computed anoma- 
lous wind, and U is the mean wind as a function of latitude. 
H' is the pycnocline anomaly computed by the model and H is 
the sloping pycnocline previously discussed. 

The key terms of equation (23) are U' and H'. Isolating 
the effect of a positive zonal wind anomaly (Fig4.6) super- 
imposed upon a negative mean zonal wind (0<U'<1), where ui = 
- 5m/sec, it can be seen that the value of (U'/U) is less 
than zero. BENZ I (ie. QO dominating all anomalies in the 
brackets), this produces a positive value of (WT) and, 
according to (12), a warming of the sea-surface temperature 
anomaly. This results because for this particular case (ie. 
Iu < 0) the total wind speed is reduced, which implies 
weaker mixing than normal. If the sign of U' and U were the 
same, then the total wind speed would be larger than normal, 
there would be greater mixing, and hence cooling of the SST 
anomaly. 

The same investigation can be done for the. physical 
effect of H'. H' is computed by the model and can be posi- 
tive or negative. In the case where H' is larger than zero, 
the total pycnocline depth is deeper than the normal B 
DH is greater than zero, and the overall sign of (WT), is 
greater than zero (see equation (23)). The interpretation 
of this result is that when the pycnocline depth is greater 
than normal, the normal wind effects cannot reach the base 
of the layer and cause the normal amount of turbulent mixing 
of cold water from below. The result is that the layer will 


be cooled less than normal.In other words, it will warm. If 
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H' is less than zero, the pycnocline is shallower than 
normal H.Then the layer will be more effected by the wind 
than normal,and therefore turbulence will cause more than 
normal entrainment, and hence the layer will experience 
anomalous cooling. It is these effects of vertical turbu- 
lent mixing which are to be investigated along with zonal 


advection (equation 12). 
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III. RESULTS AND OBSERVATIONS 


A. EFFECTS OF HORIZONTAL ADVECTION ONLY (CASE ONE) 


The control run for this thesis is the horizontal advec- 
tion ‘case  (Rennick, 1985) along with the refinements 
described in Chapter II. This means that the entrainment 
mixing term in (12) and the We' term in (13) were set to 
zero. The ocean model was integrated for a total of five 
hundred days, forced by the Blundell and Gill wind patch 
(5), with the following results. 

The pycnocline anomaly will be analyzed first due to its 
impact on the currents. Figures 4.7-4.11 show a sequence of 


"snap shots" at fifty day intervals of the current and 


pycnocline depth anomalies. The length of the arrows 
denotes the current speed in meters per second. Contour 
intervals and length scales are defined for each panel. The 


Solid contours imply a positive pycnocline anomaly and 
therefore a deeper pyncocline than normal. The dashed 
contours are negative pycnocline anomalies which indicate a 
shallower layer than normal. Figure 4.7 shows an classical 
example of a propogating Kelvin wave in that the maximum 
zonal current and maximum pycnocline anomaly locations coin- 
cide. Away from the equator,the flow is parallel to the 
isobars and is in geostrophic balance with the north/south 
pressure gradient. The movement of the pycnocline anomaly 
contours to the east is very rapid and corresponds very well 
with the propogation speed of the wind patch (approximately 
.32m/sec). 

Figure (4.8) shows the Kelvin wave spreading north and 
south along the eastern boundary. At this point, a reflec- 


muon back towards the west is evident. This continues to 
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intensify in Figure (4.9). The suface layer in the western 
ocean iS continuing to shallow, and at day three hundred 
(Figure 4.9) the shallowing has reached as far east as tage 


west longitude. 


Figure (4.10) shows the formation of an easterly jet 
along the equator. The hint of a Rossby wave to the north 
and south of the jet can be intimated. Since the Rossby 


wave moves so much more slowly than the Kelvin wave, more 
time is required to infer its actual motion. By day four 
hundred the formation of anticyclones (highs) to the north 
and south of the jet can be seen along with an intensifica- 
tion of the easterly jet caused by the Rossby waves. 

The Rossby wave motion becomes more obvious in Figure 
4.11. The westward movement of the anticyclones corresponds 
well with the Rossby wave speed. At day four hundred fifty 
(upper panel) the first indications that the pycnocline is 
Shallowing at the eastern boundary become evident. This is 
due to the the easterly jet and the cessation of the anoma- 
lous winds (see Figure 4.6). In order to better see the 
full structure of the Rossby waves, the model boundaries 
need to be more than 15° of latitude from the equator. The 
upper boundaries are causing interference, but a maximum in 
the anticyclones to the west in both panels still allows 
developement of the Rossby wave. The Rossby wave causes the 
easterly Jjet, the shallowing of the pycnocline, and the 
cooling of the sea surface from the east (see Figures 
4.15-4.16). 

Figures 4.12-4.16 are the results for SST anomaly. A 
sequence of figures at fifty day intervals for the five 
hundred days is presented. Solid contours indicate positive 
SSt anomalies and dashed contours negative SST anomalies. 
Contour intervals are specified above each panel and are in 


units of degrees centigrade. 
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Figure 4.12 demonstrates that the SST anomaly also 
propogates at the same speed as a Kelvin wave ( 2-3 m/sec.). 
By day one hundred, almost the entire basin along the 
equator has experienced some warming due to the anomalous 
advection of warm water from the west by the anomalous 
zonal currents of the Kelvin wave. An approximation in the 
Pacific is that a Kelvin wave will cross the basin in 
approximately one hundred days. This iS Supported by the 
lower panel in Figure 4.7, where it can be seen that the 
deepening of the pycnocline is already beginning at one 
hundred days. The values of the SST anomaly are very small 
and would be an observationalist's nightmare, but they are 
not random. 


Interaction with the eastern boundary is evident in 


Figure 4.13 (lower panel). The flattening of the contours 
at day two hundred implies that the Kelvin wave is 
reflecting from the boundary. In Figure 4.14 (lower panel) 


the Kelvin wave is moving north and south along the 
boundary. It should also be noted that temperatures peak 
near day two hundred fifty, and then they begin to decrease. 

The cooling from the eastern boundary continues in 
Breare 4.15, and at day four hundred the effect of the 
westerly jet caused by the Rossby waves is clearly evident. 
Figure 4.16 continues to show the effect of the easterly jet 
as cooling extends to 145° west longitude on day four 
hundred fifty. At this time the easterly wind patch ceases 
(see Figure 4.6), and cooling extends to 150° east longitude 
by day five hundred (lower panel). 

The zonal movement along the equator can better be seen 
using a time-longitude plot; therefore, Figures 4.17 and 
4.18 have been included. Solid contours are positive anoma- 
lies, and dashed contours are negative anomalies. 

The Kelvin wave front is well demonstrated in the pycno- 


cline depth anomaly at the eastern boundary by the rapid 
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increase followed by a rapid decrease (Figure 4.17). The 
decrease (dH'/dt < 0) after day three hundred is caused by 
the north/south propogation of the Kelvin wave after it has 
reached the boundary. The zonal speed of the disturbance 
may be inferred from the slope of the H'- 0 line and is 
approximately .32 m/sec, which is equivalent to the speed of 
the wind patch. 

Basin wide warming is the trend shown in Figure 4.18. 
After day four hundred negative anomalies occur due to the 
easterly jet previously discussed. The cutoff of the winds 


at four hundred fifty days is also a factor. 


B. EFFECTS OF MIXING ONLY (CASE TWO) 


For case two, the zonal advection term in the SST 
anomaly equation (12) was excluded, and the entrainment 
mixing term was included. Also, the We' term was added to 
equation (13). Similar runs were made as in case one, but 
only differences in results will be discussed here. 

The SST anomaly patterns in case two were very different 
from those in case one. The amplitude of the anomaly is very 
small but well organized. Figures 4.19-4.23 display the 
Same time "snap shots" as those done in case one Figures 
4.12-4.16, in the results. 

The SST anomaly does not propagate eastward with the 
Kelvin wave speed of 2-3 m/sec as in case one. Instead the 
progression of the SST anomaly maximum is a wind-forced 
response with a speed of approximately .32 m/sec. The meri- 
dional extent of the SST anomaly field is also much larger 
than that of case one, even though its magnitude is an order 
or two smaller. The peak SST anomaly is attained at day 
three hundred (Figure 4.21). The wind patch and its associ- 
ated SST anomaly field continues to move eastward (Figures 
4.22-4.23) until day four hundred fifty, when the wind patch 
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ceases and cooling commences at the eastern boundary. No 
negative SST anomalies ever develop in this case. The 
difference between case one (Figure 4.16) and case two 
(figure 4.23) is striking. 

The major feature in the SST anomaly field is the phase 
difference it has with the height field. The height field 
is traveling at the same speed as the SST anomaly field. 
The speed of the SST anomaly was computed by measuring the 
slope of a line through the maximums of the SST anomaly 
field in Figure 4.24. 

Rennick and Haney (1985) computed that the phase differ- 
ence should be 90° for a free wave, but in the forced case 
it is less than that. Comparing the location of the height 
field maximum (Figure 4.7, lower panel) with that of the SST 
anomaly maximum (Figure 4.12, lower panel),a phase differ- 
ence of approximately 50° was obtained. 

Although the SST anomalies in case one and two are very 
different, the pycnocline depth anomaly in the two cases are 
virtually identical. This is clearly seen in comparing 
Figure 4.17 with Figure 4.25. An important conclusion 
concerning the height field can be made from this result 
because the only difference between case one and two was due 
to the entrainment mixing parameter We' . Since no signifi- 
cant difference was evident in the two pycnocline anomaly 
fields it can be concluded that We' is negligible. Thus, 
the evolution of the pycnocline depth anomaly field is 
controlled by the dynamics of a non-entraining layer. Since 
the SST anomaly produced by anamolous entrainment mixing is 
an order or two of magnitude smaller than the advectively 
produced SST anomalies, the total SST anomaly from equation 
(12) will be almost identical to that determined by advec- 
tion alone. This conclusion would be final, if not for the 
different SST anomaly patterns produced by the two 


processes. 
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C. HORIZONTAL ADVECTION PLUS MIXING (CASE THREE) 


From the two previous cases it would be expected that 


for the linear model being used, a summation of the two 
fields will not differ much from case one. The final ques- 
tion then to ask is, "could the phase difference in case two 


have significant effect in a coupled air ocean model?".Even 
though the SST anomaly in case two is much smaller, a 
coupled atmosphere-ocean version could experience a factor 
of one hundred in the growth rate. Therefore, a one percent 
effect as seen in the uncoupled version could be magnified 
in a couple system and have an effect. This will be exam- 


ined in the following cases. 


D. .COUPLED MODE WITH HORIZONTAL ADVECTION OF SST (CASE 
FOUR) 


In case four an examination of horizontal advection in a 
coupled model is made. The question concerns whether 
growth will occur. The initial conditions for this CASCOS 
defined by the first one hundred days of case one (Figures 
4.7 and 4.12). Thus, the model is first run in a prescribed 
atmosphere mode for one hundred days, and the ocean is 
allowed to respond to the atmospheric forcing (Figures 
4.l-4.6). At day one hundred the atmosphere is released and 
allowed to respond to feedback provided by the ocean. The 
feedback of atmospheric heating is parameterized solely in 
terms of the SST anomaly predicted by the ocean model. 

Figure 4.26 demonstrates that the anomaly initially 
propogates with the wind patch, but that after day one 
hundred it rapidly increases in speed. The rapid increase 
is due to the Kelvin wave which is produced by the initial 
cut-off of the atmosphere forcing. The layer depth 
increases until approximately one hundred and fifty days 


when the Kelvin wave reaches the boundary. At this point in 
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time the entire basin begins to relax. The lack of growth 
is caused by the lack of a positive feedback. The SST 
anomaly which produces the convective heating in the atmos- 
phere simply does not develop a phase relation with the wind 
which is needed to produce a positive feedback. This result 
is consistent with the analytic findings in Rennick and 
Haney (1985) and the numerical study of Rennick (1985). 
Therefore, the conclusion can only be that horizontal advec- 
tion of SST with dynamic forcing in itself cannot cause the 
growth exhibited in ENSO. 

Figure 4.27 shows the same initial growth of SST anomaly 
as in 4.18 of case one for one hundred days. Supporting the 
findings in Figure 4.26 there is eastward propagation of the 
SST anomaly followed by relaxation of the entire basin. As 
previously stated, no positive feedback develops,and there- 


fore there is no growth. 


E. COUPLED MODE WITH TURBULENT MIXING ONLY (CASE FIVE) 


The previous case supports the hypothesis that no growth 
will occur when SST is governed by horizontal advection 
alone. Here an investigation is made concerning whether the 
different phase of the SST anomaly produced by anomalous 
mixing can tip the balance and cause growth. The model is 
initialized by repeating the first one hundred days of Case 
Three. After one hundred days the atmosphere is allowed to 
interact with the ocean, and the results are displayed in 
Figures 4.28-4.29. 

Comparing Figure 4.28 with 4.17, a large difference in 
propagation of the pycnocline anomaly is seen. Figure 4.17 
is the same for Cases One and Two and showed that the propo- 
gation speed is equal to that of the wind patch speed. This 
can be seen for the first one hundred days in Figure 4.28, 


but after one hundred days a rapid change in phase speed 
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occurs and is equivalent to the Kelvin wave speed. Thus, at 
day one hundred, a pulse is generated by the cut-off of the 
prescribed winds and it traverses the basin at the Kelvin 
wave speed (2-3m/sec). At the eastern boundary, the depth 
increases until the Kelvin wave reaches the boundary, and 
then the anomaly begins to decay. There appears to be no 
feedback or growth of any kind after one hundred and sixty 
days. The coupled system simply spins-down. 

Figure 4.29 is the SST anomaly at the equator for five 
hundred days. The difference of a factor of five in contour 
intervals also supports the results shown in Figure 4.28. 
There is no feedback or growth after one hundred days. The 
majer result is that the phase difference of the mixing case 
does not cause any growth by itself, and if added to the 
horizontal advection case, still no growth would be 


experienced. 
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IV. CONCLUSIONS 


The analysis of the anomalous entrainment mixing shows 
that in a linear ocean model and in a coupled ocean- 
atmosphere model, entrainment can be ignored as a growth 
mechanism for ENSO. The phase difference in the SST anoma- 
lies produced was not able to overcome the large difference 
between the horizontal advection case and the mixing case. 

The investigation of other terms in equation (11) must 
be investigated in order to understand fully the SST 
response to forcing fron the atmosphere in equatorial 
regions. In this analysis all that was considered was the 
anomalous entrainment mixing at the base of the layer (CWT) 
/H) and the effects of the wind on this mixing process. 

Solar radiation fluctuations have a direct effect on the 
SST anomaly as well as an effect on the turbulent mixing 
(S'(F,H) in equation 20) which has been ignored in this 
study. These effects could in turn cause anomalous heating 
or cooling of the atmosphere and have a different phase 
which may cause the growth associated with the El Nino 
event. 

The anomalous surface heat flux ( (WT), /H) would logic- 
ally be the next term to be investigated. This term is 
difficult to model because it is affected by cloud cover and 
Changes in humidity which are difficult to parameterize, 
especially in a coupled model. Therefore, a number of 
ocean-only case studies should be made, using statistical 
data from the Weare and Strub (1980) Atlas for averaged 


anomalous heat flux fields. 
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Figure 4.2 Anamolous winds for September,October, 
November 1982, observed(upper) and computer generated 
(lower); contours 5m/sec( upper) and .2m/sec(lower). 
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Figure 4.3 Anamolous winds for December,January, 
February 1982-1983, observed(upper) and computer generated 
(lower); contours 5m/sec(upper) and .2m/sec(lower). 
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a 4.4 Anamolous winds for March, April, May 1983, 
obs 


erved(upper) and computer generated (lower); 
contours 5m/sec(upper) and .2m/sec(lower). 


34 


4) » E. / Pa a ooo NAA * Y ce Ma 
ae a ` ee et qu e so r , . . a b a So - . à 
g 


— Ng S ——.— — W © 7 SS A = E Ae 


c= 9 





= ` ` ` . 





: DAE ESA 
me 0 o 


- = E - - - cams o omame PRO rn 5 
y — ~ e * ^e & £5 
o m = o = O e 


Figure 4.5 AM On winds for June, July, August 1983, 
observed( (up Hn and computer generate Molina 
contour: E RENO) and .2m/sec(lower). 
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Figure 4.6 Case one, zonal wind anomaly (m/sec) 
time-longitude plot. 
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Figure 4.8 Case one, depth (m) and current S nEs 
(m/sec) for entire basin (day 150 and 
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Figure 4.9 Case one, depth (m) and current anomalies 
(m/sec) for entire basin (day 250 and 300). 
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CURRENT AND PYCNOCLINE ANOMALIES AT DAY 350 
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Figure 4.10 Case one, depth (m) and current anomalies 
(m/sec) for entire basin (day 350 and 400). 
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CURRENT AND PYCNOCUNE ANOMAUES AT DAY 450 
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Figure 4.11 Case one, depth (m) and_current E 
m/sec) for entire basin (day 450 and 50 
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SST ANOMALY AT DAY 50 
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Figure 4.12 Case one, SST P (deg C) tor 
entire basin (day 50 and 100). 
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Figure 4.13 | Case one, SST anomalies (deg CHA EGO 
entire basin (day 150 and 20058 
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SST ANOMALY AT DAY 250 
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Figure 4.14 Case one, SST anomalies (deg E)W Tor 


entire basin (day 250 and 300 
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Figure 4.15 | Case one 
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Figure 4.16 | Case one, SST anomalies (deg Cie tor 
entire basin (day 450 and 500). 
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Figure 4.17 Case one, pycnocline depth anomaly (m) 
time- longitude plot. 
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Figure 4.18 Case one, SST anomal d E 
time- longitude plot. 7 E) 
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Figure 4.19 „Case two, SST anomalies (dee C) for 
entire basin (day 50 and 10098 
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Figure 4.20 | Case two, SST anomalies (deg C) For 
entire basin (day 150 and 200 
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Figure 4.21 ` Case two, SST anomalies (deg C) fon 
entire basin (day 250 and 300 
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Figure 4.22 Case two, SST anomalies ¡des (NEL GR 
entire basin (day 350 and 400 
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Figure 4.23 Case two, SST anomalies (deg C) fon 
entire basin (day 450 and 500 
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Figure 4.24 Case two, SST anomaly (deg C) 
time-longitude plot. 
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Figure 4.25 Case two, pycnocline depth anomaly (m) 
time- longitude plot. 
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Figure 4.26 Case four, pycnocline DD anomaly (m) 
time- Ron ELE plot 
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Figure 4.27 Case four, SST anomaly (deg C) 
time-longitude plot. 
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Figure 4.28 Case five, pycnocline depth anomaly (m) 


time longitude plot 
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